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Introduction 


The  reflection  of  a  planar  shock  wave 
over  rough  surfaces. 


Abstract 


The  effect  of  surface  roughness  on  the  transition  from  regular  (RR)  to  Mach 
reflection  (MR)  over  straight  wedges  in  pseudo-steady  flows  was  investigated  both 
experimentally  and  analytically. 

Two  models  for  predicting  the  RR  *  MR  transition  in  the  (M ,8  ) -plane  were 

s  w 

developed  (M  is  the  incident  shock  wave  Mach  number  and  6  is  the  reflecting  wedge 
s  w 

angle).  Their  validity  was  checked  against  experimental  results.  Since  tne 
experimental  results  are  limited  to  tlie  ranges  1  <  Mg  £  2  and  3.5  <  M  £  4  and  surface 
roughness  heights  of  0  <  e  £  0.2  cm,  the  proposed  models  are  applicable  bo  these 
ranges  only.  In  the  first  model  (the  pressure  loss  model' ,  the  transition '  from  PJR 
to  MR  is  assumed  to  be  relate':  to  tlx?  pressure  loss  due  to  the  increased  friction 
oetveen  the  flow  and  the  rough  surface.  In  tix2  second  model  (tnc  boundary  layer 
displacement  thickness  model)  ,  trie  HR  t  MR  transition  is  related  to  the  boundary  layer 
displacement  thickness  which  in  turn  depends  on  the  surface  roughness. 
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It  was  snown  that  the  boundary  layer  displacement  tnickness  6*  could  be 

6LV 

characterised  by  the  parameter  G  =  — (where  <5^  is  the  laminar  sublayer 
thickness,  V*  is  the  shear  velocity  and  v  is  the  kinematic  viscosity) .  This 
parameter  was  found  to  be  constant  for  a  given  size  of  surface  roughnffis- e  i.e.,  G  was 
found  to  be  independent  of  the  incident  shock  wave  Mach  number  -  Mg  .  Owing  to  this 
fact  G  was  termed  the  "roughness  character izer." 
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Introduction 

When  a  planar  incident  normal  shock  wave  encounters  a  sharp  corner  in  a  shock 

tube,  four  different  types  of  reflection  can  occur  depending  on  the  incident  shock 

wave  Mach  number  -M  arc  the  wedge  angle  -  6,..  They  are : 

s  w 

1)  regular  reflection  -  HR  (figure  la) , 

2)  single-Mach  reflection  -  SMF  (figure  lb) , 

3)  complex-Mach  reflection  -  G4R  (figure  lc) ,  and 

4)  double-Mach  reflection  -  DM.R  (figure  Id) . 

The  latter  three  types,  i.e.,  SMR,  G-1R  and  DMR  are  usually  termed  as  .Mach  reflections 
-  MR 

The  RR  *  MR  transition  criterion  in  pseudo-steady  flows  (the  "detachment 
criterion")  was  first  introduced  by  von  Neumann  (1963) .  In  developing  this  criterion 
von  Neumann  assumed  that: 

1)  the  flow  is  ideal,  i.e.,.  w  =  0  and  k  *  0  ( v»  is  the  dynamic  viscosity  and  k  is 
the  heat  conductivity) , 

2)  the  flow  is  two-dimensional,  and 

3)  the  flow  is  self-similar  and  hence  pseudo-steady. 

As  shown  by  Ben-Dor  &  Glass  (1979),  the  detachment  criterion  can  be  expressed  as 
follows: 

61  +  e2m  =  0  ll> 


where  6^  is  the  flow  deflection  through  the  incident  shock  wave  -  i,  and  6^  is  the 
maxima”,  possible  flow  deflection  through  the  reflected  shock  wave  -  r.  This 
formulation  is  uased  on  the  fact  that  when  the  frame  of  reference  is  attached  to  the 
reflection  point  -G  (see  figure  2)  the  flow  ahead  of  the  incident  shock  wave  -  i, 
whicti  now  moves  parallel  to  the  wedge  surface  at  supersonic  velocity  v  =  u  /oos9 


w 


(  u  is  tne  incident  snock  wave  velocity) ,  is  deflected  towards  the  wedge  by  on 


angle  6^  .  Thus  tne  supersonic  flow  behind  the  incident  shock  wave  -  i  is  deflected 


towards  the  solid  surface.  In-order  to  negotiate  the  solid  surface  it  must  experience  a 

reeeflection  while  crossing  tne  reflected  snock  wave  r  to  become  again  parallel  to 

the  wedge  surface.  When  the  angle  0  exceeds  the  maxi mux  deflection  angle  0  ,  the 

2  2m 

regular  reflection  becomes  theoretically  impossible.  Consequently  equation  (1) 
represents  the  "detachment"  criterion  of  von  Neumann. 

However,  the  experimental  findings  of  von  Neumann  (1943),  Smith  (1945),  Taub 
(1947),  Bleakney  &  Taub  (1949),  Fletcher  (1950),  VJhite  (1951),  Law  &  Class  (1971), 
Henderson  &  Lozzi  (1975),  Ben-Dor  (1978),  Henderson  &  Woolington  (1983)  and  Hornung  & 
Taylor  (1982)  indicates  tnat  the  RR  exists  beyond  the  predicted  theoretical  limit  of 

the  "detachment"  criterion. 

The  transition  line  predicted  by  the  "detachment"  criterion  (equation  1 ) ,  as 
veil  as  our  experimentally  obtained  RR  t  MR  transition  points  are  shewn  in 
figure  3.  The  persistance  of  the  RR  beyond  the  limit  predicted  by  the 
detachment  criterion  is  clearly  seen.  Since  the  von  Neumann  paradox  was 
first  noticed,  many  scientists  (Henderson  &  Lozzi,  (1975) ,  Ben-Dor  &  Glass  (1979) , 
Henderson  &  WooLmington  (1983),  Hornung  &  Taylor  (1982)]  directed  their  efforts  at 
resolving  this  paradox.  Probably  the  most  promising  approach  was  tne  one  suggested  by 
Hornung  f«  Taylor  (1982).  They  argued  that  the  re  .son  for  the  existence  of  the  von 
Neumann  paradox  is  tne  fact  that  the  transition  line  was  derived  by  solving  the 
inviscid  flow  conserv.it ion  equations  while  the  actual  flow  was  viscous.  Consequently, 
they  concluded  that  by  accounting  for  viscous  effects  the  von  Neumann  paradox  could  be 
resolved.  Hornung  f.  Taylor  (1982)  accounted  for  the  viscous  --ffects  by  applying  tne 
ooundary  layer  displacement  thickness  concept.  Tne  obtain.:  ;  results  justil. Lv,  their 
approach. 

Takayama,  Gotoh  and  Ben-Dor  (1981)  presentee  ox peri. to  nt.il  results  regarding  hi  £ 

V.P  transition  over  rough  surfaces.  Their  results  are  reproduce:  in  figure  4.  Tne 
difference  between  the  actual  transition  anti  the  values  predicted  by  the  detachment 
criterion  (solid  line)  is  quite  clear.  Tins  difference  increases  as  the  roughness 
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neight  increases.  For  example,  at  Mg  =  4  and  e  =  0.2  cm,  tne  actual  wedge  angle  at 
which  t'ne  RP  *  MR  transition  occurs  is  about  20°  lower  than  that  predicted  by  the 
"detachment"  criterion. 

The  fact  that  the  surface  rougnness  has  such  a  meaningful  influence  on  the 
transition  wedge  angle  on  one  hand,  and  the  fact  that  shock  reflections,  in  nature, 
occur  over  rough  surfaces,  has  led  to  the  recognition  that  understanding  the 
reflection  process  over  rough  surfaces  is  of  great  importance.  Consequently,  models 
capable  of  predicting  the  RR  -e  MR  transition  over  rough  surfaces  were  sought.  This 
was  further  stressed  in  the  first  concluding  remark  of  the  3rd  Mach  Reflection 
Symposium  (held  in  Melbourne,  Austrlia,  in  August  1983)  whicn  said:  "It  appears  to  be 
the  ooundary  layer  which  causes  the  delay  of  the  transition  from  regular  to  Mach 
r<.  flection  of  pseudo-stationary  shocks  beyond  the  detachment  or  maximum  deflection 
point  for  reflection  from  a  plane  surface.  This  effect  is  more  pronounced  in  the  case 
of  rough  surfaces." 

Based  on  the  foregoing  discussion  the  RR  *  MR  transition  process  was  investigated 
botn  experimentally  and  analytically.  Since  Hornung  &  Taylor  (1982)  were  able  to 
explain  the  von  Neumann  paradox  by  including  the  viscous  effect  in  their  solution  of 
the  flow  field,  it  was  decided  to  base  the  present  models  for  predicting  the  RR  t  MR 
transition  over  rough  surfaces  on  "viscous  mechanisms." 

Subsequently,  two  different  approacnes  tnrough  which  the  effect  of  viscosity  can 
be  accounted  for  will  oe  proposed.  In  the  first,  tiie  pressure  loss  due  to  friction  is 
integrated  into  the  invisciu  conservation  equations,  while  in  t'ne  second  the  wedge 
geometry  is  modified  using  the  boundary  layer  displacement  thickness  concept. 

1  xper imental  Results 

Tne  experimental  results  re[>orte<.:  here  were  obtained  using  tin-  40  x  80  mm  shock 
tube  of  the  IHSM  Institute  of  High  5;peed  Mechanics,  Tohoku  University,  Sendai,  Japan. 
An  illustration  of  the  "saw-tooth"  roughness  imposed  on  the  tested  wedge  models  is 
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given  in  figure  5a.  For  the  first  set  of  experiments  four  models  (each  4  cm  long) 

with  roughness  of  c  =  0.01,  0.02,  0.08  and  0.2  cm  were  prepared.  The  effective  wedge 

angle  of  the  models  could  be  adjusted  frctn  outside  the  shock  tube  test  section  as  shown 

in  figure  5b.  The  experimental  study  covered  the  range  of  1  <  M  £  4. 

s 

The  shock  wave  velocity  was  measured  using  Kistler  606  pressure  tranducers,  an 

Iwatsu  UC  7641  digital  counter  and  an  Iwatsu  DM  901  recorder.  The  pressure  tranducers 

were  mounted  120  mm  apart  just  ahead  of  the  test  section.  The  attenuation  of  the 

incident  shock  wave  was  found  to  be  negligibly  small. 

The  test  gas  was  pure  nitrogen.  The  initial  pressures  were  in  the  range  of  3  £  Pq  £ 

150  torr.  For  the  investigated  range  of  Mach  numbers  (1  <  Mg  £  4)  the  nitrogen  gas 

can  be  assumed  to  behave  as  a  perfect  gas.  Therefore,  the  initial  pressure  has  no 

effect  on  the  reflection  phenomenon.  For  all  the  experiments  the  initial  temperature 

was  about  T  =  300K. 
o 

An  optical  diagnostic  (shadowgraph)  was  used  for  studying  the  pseudo-steady 
reflection  phenomenon.  A  pulsed  ruby  laser  (6943  A)  was  used  as  a  light  source. 

Typical  shadowgraphs  are  shown  in  figures  6a  to  6d.  Mach  reflection  over  a  wedge 
having  roughness  height  of  e  =  0.02  cm  is  shown  in  figure  6a.  Reqular  reflection 
over  the  same  wedoe  is  shorn  in  figure  6b.  Ficrrues  6c*  and  6d  show  a.  MR.  and  a  HR 
over  a  wedge  having  roughness  height  of  e  =  0.2  an,  respectively. 

The  W  *  Mfc  transition  angle  was  determined  in  the  following  way.  For  a  given 

value  of  tiie  incident  shock  wave  Mach  number  -  M  the  wedge  angle  e  was  gradually 

s  w 

increased  until  a  RR  was  obtained.  Tnen,  the  value  of  the  triple  point  trajectory 
angle  x  (deduced  from  the  snadowgraphs)  was  plotted  in  the  (x»9w  ) -plane  as  shown  in 
figure  7.  The  experimental  results  were  then  extrapolated  to  x  =  0.  The  wedge  angle 

at  tins  point  was  cnosen  as  the  transition  wedge  angle  -  e  .  This  procedure 

v? ,  tr 

yields  0^  to  an  accuracy  of  ±1  “ .  The  above-mentioned  procedure  was  repeated  for 

Mg  =  1.04,  1.12,  1.21,  1.44,  1.96,  3.58,  3.77  and  3.89. 

Tlie  experimental  results  obtained  for  the  four  different  models,  as  well  as  the 
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! 

i 

results  obtained  for  a  "smooth"  model,  are  shown  in  figure  4.  * 

J 

It  is  apparent  that  the  surface  roughness  (£)  has  a  significant  influence  on  the  j 

I 

RR  t  MR  transition  wedge  angle  for  any  given  incident  shock  wave  Mach  number.  The  [ 

I 

greater  e  is,  the  smaller  is  the  transition  wedge  angle  6w  tr»  It  is  of  interest  to  i 

note  that  the  measured  transition  wedge  angles  over  the  smooth  surface  model  does  not 
agree  with  the  "detachment"  criterion  (these  results  confirm  the  von  Neumann  paradox) . 

Had  the  surface  been  perfecty  smooth,  i.e.,  £  =  0,  the  experimental  results  should 

i 

have  agreed  with  the  "detachment"  criterion.  However,  since  no  surface  is  perfectly 
smooth,  the  von  Neumann  paradox  exists.  In  the  following  we  will  refer  to  a  "smooth" 
surface  as  hydraulically-smooth.  Furthermore'',  it  will"  be  proven  analytically  that 
a  hydraulically  snooth  .surface  is  a  surface  for  which  e  g  0.00517  an. 

Analysis 

As  mentioned  previously  the  aim  of  the  present  study  is  to  develop  models  capable 
of  predicting  the  RR  t  MR  transition  over  rough  surfaces,  i.e.,  models  which  will 
shift  the  "detachment"  transition  line  (solid  line  in  figure  4)  towards  the 
experimental  results  for  any  given  roughness  height  of  -  e  . 

1 

i 

General  Assumptions  ! 

j 

1)  The  flow  is  two-dimensional. 

2)  The  flow  field  can  be  described  by  the  mass,  momentum  and  energy 
conserrvation  equations  of  inviscid  flows  with  additional  terms  accounting  for 
viscous  effects. 

3)  The  flow  is  self-similar,  and  hence,  can  be  made  psoudo-stoady  by  applying 
tiie  well-known  Galilian  transformation. 

4)  Real  gas  effects  can  be  ignored,  consequently  the  gas  is  assumed  to  oliey  the 
equation  of  state  of  a  perfect  gas. (i.e.,  P  =  pKT) . 
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The  Conservation  equations 

The  wave  configuration  of  a  PJv  is  shown  schematically  in  figure  2.  linploying  the 
foregoing  assumptions,  and  attaching  the  frame-  of  reference  to  the  reflection  point  6 
(in  such  a  way  that  it  moves  along  the  wedge  surface  together  witn  point  G)  results  in 
the  following  nine  equations  which  completely  describe  the  flow  field  in  the  vicinity 
of  the  reflection  point  -  G. 
across  the  incident  shock  wave  i: 


continui ty 


°ovosin*o  =  P1v1s^Vel) 


tangent ional  momentum 


potai^o  =  P^VV 


normal  momentum 

po  +  povo2sin\  =  pi(1~V  + 


energy 


ho  +  \  Vo2sijl\  =  hl  +  I  v^sii^-e,) 


across  the  reflected  snock  wave  t: 


(2) 


(3) 


(4) 


(5) 


continuity 

P1v1sin<j.1  =  p2v2sin(4>1-02) 


tangent ional  momentum 

P1tan^1  =  p2tan(<j>1-e2) 


(6) 


(7) 
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normal  momentum 

P1  d-Fp)  +  p1v12sin24>1  =  P2  (1— Fp)  +  p2v22sin2  (4>1“62)  (8) 


energy 

12  2  12  2 
h1  +  2  V1  sin  *1  =  ^2  +  2  V2  Sin  ^1_02^ 


(9) 


and  the  boundary  condition  for  the  RR  t  MR  transition: 

•,*%  =  “  (10) 

Setting  =  0  and  a=  0  results  in  the  well-known  set  of  equations 

describing  a  regular  reflection  in  an  inviscid  flow.  The  solution  of  tnese  equations, 

which  yields  the  "detaciment"  transition  line  (shown  in  figures  3  and  4),  can  be  found 

in  Ben-Dor's  (1978)  work.  Note  that  v  =  u  /cose  ,  6  =90-0  • 

o  s  w  yo  w 

and  Tq  are  known.  The  enthalpy  h  is  related  to  the  tempertaure  through  h  =  CpT 
and  the  density  P  can  be  calculated  from  the  equation  of  state  for  a  perfect  gas, 
i.e.,  p  =  P/RT  .  The  definitions  of  <f>Q,  ,  6^  and  fire  shown  in  figure  2. 

The  Roughness  Effect  on  the  Turoulent  Liow 

Usually,  surfaces  over  which  fluids  flow  cannot  be  considered  as  perfectly 
smooth,  hxperiments  have  indicated  tnat  the  roughness  height  -e  significantly 
contributes  to  tne  friction  coefficient  -  f.  Owing  to  tnis  fact  it  is  necessary  to 
fiovtlop  models  whicn  will  relate,  the  friction  coefficient  f  to  the  roughness 
height  -  e  .  'mis  is  done  in  tne  following. 

Llasius  [behl ichting  (1902) )  snowed  experimentally  that  for  hywraulically  sr.ooth 
pipes  tne  friction  coefficient  f  do[iends  solely  on  the  Reynolds  number  (Re),  liopf 
(1923)  found  experimentally  tnat  for  very  rougn  pities  the  friction  coefficient  f 
depends  solely  on  the  relative  roughness  e/R  (R  is  the  pipe  radius) .  Furthermore,  he 
si  lowed  that  for  less  rough  pipes  f  depends  both  on  Re  and  g/R.  These  experimental 
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firriings  suggested  that  the  considered  niienanenon  is  controlled  by  tin-  thickn-  of 
too  laminar  sublayer  -  6  .  Their  conclusion  was  t.int  trie  dominant  factor  in 
determining  the  value  of  tne  friction  coefficient  f  is  e/6  .  alien  c  <  5  (i.t.,  tne 

surface  roughness  is  completely  immersed  inside  tne  laminar  sublayer)  the-  roughness 
has  practically  no  effect  on  £,  and  the  surface  can  be  considered  as  "hydraulically 

smooth."  On  the  other  hand,  for  e  >>  6  the  roughness  has  a  dominant  effect  on  f.  In 

L 

such  a  case  the  friction  mechanism  is  known  as  "wave  drag"  [L-hames  (1982)].  It  was 


further  shown  by  Schlichting  (1962)  that  for  any  given  size  of  roughness, 


V*6L 


const  (v*  is  the  shear  velocity  and  v  is  the  kinematic  viscosity).  This  finding  led 
Sclicnting  to  derive  the  following  relation: 


ev* 

v 


(ID 

It  is  clear  from  this  relation  that  -  can  replace  ■£—  and  therefore  lie 

v  6t 


ev* 


considered  as  the  dominant  factor  in  determining  the  value  of  the  friction 
coefficient. 

Basea  on  the  experimental  results  of  Nikuradsse  (1933)  it  is  a  common  practice  to 
divide  tne  turbulent  flow  over  rougti  surfaces  into  three  flow  zones  [Shames  (19s2) | . 


ev 


1)  hydraulically  smooth  flow:  —  <  1  or 

6L 


<  5  and  f  =  f  (Re) 


ev* 

2)  frictional  transition  flow:  1  <  ~~  <  15  or  5  <  -  <  70  and  f  =  f(Re,  e/R) 

5l  v 

eV* 

i)  cough  flow:  ~  >  15  or -  >  70  and  f  =  f (e/R) 

6L 


Oik.'  tne  turbulent  character  of  the  flow  is  defined,  it  is  important  to  ivv« ■ 
correlations  relating  the  roughness  to  the  friction  losses.  As  a  basis  for 
constructing  seen  correlations  the  following  logarithmic  velocity  distribution  law  is 


uses : 


u  1  M 
rr-  =  —  I  In 


v*y 


-  lne] 


(12) 
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where  u  is  the  time  average  veloci'y,  is  a  universal  constant,  y  is  tn-_  distance 

coordinate  mc-osured  from  the  surface  and  6  is  a  constant  depending  on  tne  roughness 
height  -  e.  The  assumptions  under  wnich  the  loyar ltmic  velocity  distr ioution  law  was 
obtained  are  given  by  Prandtl  (1925).  Detailed  derivation  of  equation  (12)  can  be 
found  in  Shames  (1982)  .  A  comparison  between  the  experimental  results  of  t.'ikuraaso 
(1933)  and  equation  (12)  indicates  that: 

1)  for  a  large  range  of  Reynolds  numbers  and  roughness  lieights,  =  0.4 


2)  for  a  hydraulically  smooth  flow 
In 6  =  -2.2 

for  the  frictional  transition  flow, 

-2.2  *  lne  <  0.8 

and  for  the  rough  flow 

eV* 

in  8  =  In - 3.4 

v 


(13) 


(14) 


(15) 


Based  on  these  values  of  and  8  expressions  for  evaluating  t'e  friction 


coefficient  -  f,  for  flat  plates,  were  derived  for  each  of  trie  thteo  flo*  /on.  s 
order  not  to  limit  the  present  results  to  any  spec'ifio  flow  zone  a  ~.ifler«.: 
analytical  approach  is  taken  in  the  present  study.  This  approuca  is  rut  lined  in  t 
following . 


J  n 


Inserting  the  boundary  condition  u  ^ 


(figure  8)  into  equation 


(12)  yields 


Ur  1  V*6t 

^  =  ^  ^  ^  ^ 


In  the  laminar  sublayer  the  velocity  changes  linear i 3 y  therefore 

Tw  .  ' 

°L  =  T  6l 


(1  ) 


(17) 


2 

where  tw  is  the  shear  stress  on  the  surface.  Inserting  the  relation  v*  =  t p 


«L  V*6I 


V*  v 


(16) 


Combining  equations  (16)  and  (18)  yields: 


1 


G  =  —  (lnG-lnB) 

v*6l 

where  G  =  — - — 


(19) 


For  =0.4  equation  (19)  reduces  to 


G  =  2.5(lnG~ln8) 


(20) 


Since  6  depends  solely  on  the  roughness  height  e,  it  is  obvious  from  equation  (20) 
that  so  Joes  G,  i.e.,  for  any  given  value  of  e,  the  value  of  G  is  independent  of  the 
flow  properties.  It  was  therefore  decided  to  call  G  the  "roughness  character izer". 


VL 

It  should  be  noted  that  the  roughness  characterizer  G  =  — —  is  sinply  the 


ocmnon  Reynolds  nunber  in  the  usual  form  applicable  to  the  inner  layer.  Hie 
laminar  sublayer  thickness  <5^  is  fundamental  in  the  evaluation  of  G.  Indeed,  G  does 
characterize  the  roughness  effects,  but  should  more  properly  be  viewed  as  the  ratio 
of  two  length  scales;  6L  and  a  viscous  length 

For  assessing  the  value  of  G  appropriate  to  hydraulically  smooth  surfaces, 
Blasius'  results  can  be  used.  Blasius  found  [Schlichting  (1982)]  that  the  friction 
coefficient  fo:  tne  considered  case  is  given  by: 


f  =  0.3164  R 


-1/4 


-acre  R 


u  D 
av 


le,D  (21) 

.  Equation  (21)  ,  which  is  known  as  the  Blasius  equation,  was 


e,D  v 

derived  for  the  following  specific  velocity  profile: 


1 


—  =  <2)7 

u  'r; 
o 


(22) 


For  a  general  velocity  profile. 


iL  •*-(*) 
uo  R 


n 


equation  (21)  can  be  generalized  to  the  following  form: 


-  2n 

f  =  C(n,G)  Re 


v  ,-N. ... . 


(23) 


In  this  relation  the  velocity  power  exponent  n  depends  solely  on  tne  Reynolds 


number  (Shanes,  1982) .  Mazor  (1984)  shewed  that  CXn,G)  has  the  following  form. 


5n+1  2  2  (n+1) 

G(n,G)  =2  n+1  [  (n+1)  (n+2)  ]n+1  G  n+1 


(25) 


Inserting  n  =  into  equation  (25)  and  recalling  that  for  this  value  of  n, 
C(n,G)  =  0.3164  (see  equation  (21)]  yields  for  a  hydraulically  smooth  surface 


G  =  12.2468 

Using  this  value  of  G  in  equation  (20)  implies  that  for  a  hydraulically  smooth  surface 

ln8  =  -2.39 


This  value  of  InB  is  about  8%  smaller  than  the  value  usually  quoted  in  the 
literature  [equation  (13)]. 

The  value  of  the  roughness  characterizer  -  G  for  a  frictional  transition  flow 
can  be  determined  from  equation  (20).  As  mentioned  earlier,  for  a  frictional 
transition  flow  -2.2  $  InB  $  0.8  [equation  (15)],  thus  by 

inserting  InB  =  0.8  into  equation  (20),  the  lower  value  of  G  for  a  frictional 
transition  flow  can  be  obtained. 

Following  this  procedure  one  finally  obtains: 

-  for  the  hydraulically  smooth  flow  G  =  12.2468 

-  for  the  frictional  transition  flow  0.5625  <  G  <  12.2468 

-  for  the  rough  flow  0  <  G  <  0.5625 

Models  for  Predicting  the  Effect  of  the  Surface  Roughness  on  the  RR  t-  MR  Transition 

As  mentioned,  the  detachment  criterion  of  von  Keumann  (1963)  Joes  not  agree  with 
experimental  results.  By  accounting  for  viscous  effects,  llornung  i.  Taylor  (1932)  were 
able  to  explain  the  discrepancy. 

For  rough  surfaces,  where  viscous  effects  are  much  more  pronounced,  the  inviscid 
model  (equations  (2)  to  (10)]  predict  the  RR  t  MR  transition  very  poorly.  In  tne 
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following,  two  models  for  predicting  the  RR  t  Hl<  transition  over  rough  surfaces  are 
developed.  The  formulation  of  these  models  is  similar  to  the  inviscid  flow  model, 
with  additional  terms  accounting  for  viscous  effects. 

The  Pressure  Loss-Model 

Viscosity  results  in  pressure  losses.  Therefore  it  is  natural  to  express  its 
effects  on  the  flow  field  by  introducing  a  pressure  loss  term  in  the  inviscid 
conservation  equations.  Since  ahead  of  the  incident  shock  wave  the  flow  is 
quiescent,  only  the  pressures  in  states  (1)  and  (2)  [see  figure  2]  need  to  be 
corrected.  Consequently,  a  pressure  loss  coefficient  Fp  ,  which  depends  on  both  the 
roughness  height  e  and  the  inverse  pressure  ratio  across  the  incident  shock  wave  -  £ 
is  introduced  into  the  conservation  equations  for  inviscid  flow.  This  pressure  loss 
coefficient  can  be  expressed  as 

AP  =  Fp(e,d  (P1“P2)  (26) 

Tnus  equations  (2)  to  (10)  with  Fp^0  and  a  =  0  represent  the  conservation 
equations  for  the  pressure  loss  model. 

Takayama,  Gotoh  &  Ben-Dor  (1981)  showed  that  the  pressure  loss  coefficient  can  be 
expressed  as  f Hows: 

Fp«c,e;»  =Fc<1-5)*'  (27) 

wnere  F^  depends  on  the  roughness  height  c  only  and  n  is  a  constant;  n  can  be  found 
from  the  experimental  results.  Their  exper imentol  data  indicate.!  that  n  =  7/9.  The 
RK  t  HR  transition  lines,  as  predicted  by  equations  (2)  to  (10)  for  Fp^O  and  a  =  0 
were  fitted  to  the  exper imentol  results  of  Takayama,  Goton  and  Ben-Dor  (1981)  by 
cnoosing  appropriate  values  for  F^  .  The  obtained  results  are  shown  in  figure  9. 

The  values  of  F^  wnich  resulted  in  the  desired  shift  of  the  RR  J  HR  transition  lines 
(to  fit  the  experimental  results)  are  listed  in  figure  9.  It  is  apparent  from  figure 
9  that  the  experimentally  obtained  transitions  for  hydraulically  smooth  surfaces  can 
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best  be  described  by  F^  =  0.1.  The  value  =  0  corresponds  exactly  to  the 
"detacnment"  transition  criterion  of  von  Neumann,  [i.e.,  a  perfectly  smooth  surface 
(  e-  0)]  . 

In  order  to  complete  this  model,  the  values  of  F^  were  correlated  to  tne 
roughness  height 

2 

F£  «  -  545.45  e  +  31.73c  -0.012;  for  0  *  e  £  0.02  cm  (28) 

and 

F£  =  -4.63e2  +  2.55e  +  0.35;  for  0.02  S  e  $  0.2  cm  (29) 

These  correlations  along  with  the  measured  values  of  ( F  -  e  )  ore  shown  in 
figure  10. 

In  summary,  the  proposed  model  should  be  used  as  follows.  For  any  given  value 
surface  roughness  c  (in  cm)  one  evaluates  using  equations  (28)  or  (29).  Then  Fp 
is  calculated  by  equation  (27) .  Subsequently  this  value  is  used  in  the  conservation 
equations  ((2)  to  (10)].  It  should  be  noted  that  for  the  pressure  loss  model  a 
[equation  (10)]  should  be  kept  equal  to  zero. 

The  Boundary  Layer  Displacement  Thickness  Model 

The  idea  to  account  for  the  viscous  effects  through  the  boundary  layer 
displacement  thickness  is  due  to  Hornung  f.  Taylor  (1982)  who  resolved  the  von  Neumann 
paradox  by  applying  tnis  concept. 

The  ixxindary  layer  displacement  thickness  -  6  *  is  defined  as  the  distance  by 
which  the  solid  boundary  would  nave  to  be  displaced  if  the  entire  flow  was  replaced  by 
a  similar  frictionless  flow  both  having  the  same  mass  flow  at  any  flow  cross-section. 

aince  the  incident  shock  wave  was  brought  to  rest  by  using  a  Galilian 
transformation,  the  boundary  layer  displacement  thickness  in  the  present  case  is 
negative  (see  figure  11) . 

In  reality  the  flow  associated  with  the  RK  is  viscous,  therefore  solving  it  as  an 
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mviscid  flow  (figure  12a)  is  inadequate.  If  one  is  to  use  the  inviscid  conservation 
equations,  the  geometrical  boundary  must  oe  corrected  in  order  to  account  for  the 
viscous  effects.  One  possible  way  is  by  displacing  the  flow  boundary  (wall)  according 
to  the  boundary  layer  displacement  thickness.  The  resulting  flow  field  is  shown  in 
figure  12b.  As  can  be  seen  the  boundary  condition  for  RR  is  no  longer  6^  +6^=0. 
The  boundary  condition  that  the  flow  obeys  behind  the  reflected  shock  wave  is 
0^+02=a  where  a^O  .  Therefore,  at  transition 


or  alternately. 


61  + 


0~  =  a 

2m 


1  + 


3  2m 


(30) 


(31) 


where  F  -  the  boundary  layer  displacement  thickness  coefficient  is  defined  as 

O 

=  a/0^  .  Consequently,  the  solution  of  equations  (2)  to  (10)  with 

Fp  =  0  and  a  ^  0  should  provide  the  correct  RR  t  MR  transition  for 

rough  surfaces  since  it  accounts  for  the  surface  roughness  via  the  displacement  thickness. 

In  the  following  the  relation  between  F.  and  e  is  developed.  It  is  clear  from 

6 

figure  4  that  for  any  given  value  of  roughness  size  e,  the  actual  transiton  points 
shift  away  from  the  "detaclment"  criterion  transition  line  to  smaller  wedoa 

angles.  Thus,  by  selecting  an  appropriate  value  for  F  the  transition  line  could  also 

6 

be  shifted  to  fit  tne  experimental  results. 

The  required  values  of  F  for  tne  given  values  of  M  aid  e  are  listed  in  table 

o'  s 

1.  Inspecting  these  values  indicates  that  in  the  range  1  <  K  <:  2  and  3.5  £  N.  <  4.0  1' 

s'  's'  6 

is  practically  independent  of  M  ,  i.e.,  F  cun  be  assumed  to  depend  soley  on  e. 

S  0 

Owing  to  this  finding,  in  the  subsequent  discussions  two  separate  ranges  of  M  are 

"  «P 

considered,  1  <  M  £  2  and  3.5  <  M  <  4.  The  values  chosen  for  F  for  the  range  1  < 

s  s  6 

j:  2  are  shown  in  figure  13.  It  is  apparent  from  figure  13  that  the  selected 
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values  of  F,  reproduce  the  experimentally  obtained  RR  t  MR  transitions  for  any 
6 

roughness  size  e  fairly  well. 

The  values  of  f.  and  a  which  reproduce  the  experimental  results  for  the  two 
6 

ranges  of  the  incident  shock  wave  Mach  number  ranges  are  given  in  table  2. 

The  correlations  between  F  and  e  ,  based  on  the  experimental  results  shown  in 

6 

figure  13,  are: 


for  1  <  M  <2 

s  - 

F  =  -600  £2  +  38.6  e  +  0.006 

6  e 

F6  =  -4.17  e2  +  3.08  e  +  0.48 


and  for  3.5  £  M  4.0 


0  <  e  £  0.02  cm  (32) 
0.02  £  c  £  0.2  cm  (33) 


F  =-9.45.45  e  +  37.13  e  -0.003  0  <  e  S  0.02  cm  (34) 

F  =  -1.39  e2  +  2.47  e  +  0.31  0.02  £  e  £  0.2  cm  (35) 

Obviously  for  a  perfectly  smooth  surface  fe  =0)  viscous  losses  can  be  ignored  and 
therefore, 


F6  =  0  (36) 

Equations  (32)  to  (35)  along  with  the  experimentally  deduced  values  of  F^  and  e 
are  plotted  in  figure  14. 

In  summary,  the  procedure  for  employing  the  boundary  layer  displacement  thickness 
model  is  as  follows:  for  a  given  value  of  roughness  height  -e  (in  cm)  and  incident 
shock  wave  Mach  number  -  the  value  of  F^  is  calculated  from  the  appropriate 
correlation  (equations  (32)  to  (36)].  Then  the  conservation  equations  ((2)  to  (1U)) 
are  solved  using  the  obtained  value  of  F^  (with  F  =  0) .  besides  its  ability  to 
reproduce  the  RR  Z  MR  transitions  for  rough  surfaces,  the  proposed  model  (the  boundary  layer 
displacement  thickness)  can  be  justified  on  physical  grounds.  It  will  be  shown  in  the 
following  tnat  for  a  given  value  of  roughness  height  -  e,  the  slope  of  the  boundary 
layer  Displacement  thickness  -  6*  at  a  specified  point  (which  will  be  termed  "the 
characteristic  distance  xf,har  ")  is  equal  to  the  angle  a  which  can  lie  obtained  from 
the  model  described  earlier  Much  is  based  on  experimental  results)  for  the  same  value 
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of  e.  (Recall  that  a  =  F„/8^  is  the  angle  between  the  flow  direction  behind  the 

o  1 

ref lecterl  shock  wave  and  the  real  reflecting  surface  (see  figure  12b)),  i.e., 


a  =  tan 


-1  d6* i 


Intuitively,  one  might  suggest  that  since  the  displacement  starts  at  the  reflection 


point  G  (figure  12b),  equation  (37)  should  be  correct  at  *char  =  0.  However,  since 


-1  d6*, 


and  since  a  flow  cannot  follow  a  90°  turn  instantaneously 


it  appears  that  equation  (37)  must  hold  for  x 


0.  Furthermore,  following 


Shirouzo  &  Glass  (1982)  and  Mireles  (1983)  it  is  preferable  to  use  the  average  slope 
of  the  displaced  boundary.  Thus,  equation  (37)  is  replaced  by: 

.  -1  6*. 

a  =  tan  — 

x  l^char  (38) 

In  order  to  proceed  with  the  present  analysis,  an  expression  for  the  boundary 
layer  displacement  thickness  -6*  of  a  compressible  viscous  flow  over  rough  flat  plates 
is  required.  An  appropriate  expression  is  developed  in  appendix  A.  The  obtained 
expression  for  6*  (see  equation  A. 17  in  appendix  A)  is: 

2(n-1)  m  3.2975n  _2n_ 

6*  _  ^  3n+1  ,V  3n+1  ,  2  „  3n+1  2, 

-  °=  <tT»  (1v,^  I*  '  '  T"’  <39> 

2  I  2  o|  m 

where  C  is  a  constant  depending  solely  on  n  (equation  A. 18),  n  is  the  velocity  profile 
exponent,  which  depends  solely  on  the  Reynolds  number  [Shames  (1932),  p.  281],  G  is 
the  roughness  character izer  [defined  by  equation  (19)],  T2,  V2  and  V2  are  the  flow  temp 
erature,  the  flow  velocity  and  the  kinematic  viscositv  behind  the  reflected  shock  wave 


respectively,  T^  is  the  average  temperature  inside  the  boundary  layer  [defined  by 
equation  (A. 8)]  and  Vq  is  the  plate  velocity  on  which  the  boundary  layer  develops. 


Inserting  equation  (39)  into  equation  (38)  yields: 


f  2 (n-1 )  3.2975n 

.  -1  J  __  3n+1  3n+1  . 

a  -  tan  1  OG  (—1  (-1 


-)3n+1  (n+1 


There  are  three  unknown  parameters  in  the  rignt  hand  side  of  equation  (40);  they 


are: 

1)  the  roughness  character izer  -  G  which  in  the  proposed  model  depends  only  on 
the  roughness  height  -  e  (this  fact  is  yet  to  be  proved) , 

2)  the  velocity  profile  exponent  -  n  which  according  to  Shames  (198?)  should 
decrease  with  increasing  Reynolds  numbers  (this  dependence  is  yet  to  be  found) , 


3)  the  characteristic  distance  -  ,  which  according  to  the  present  model 

snould  approach  zero,  i.e.,  xr,ywr  •+•  0* 

In  the  following  these  unknowns  will  be  found.  However,  in  order  to  do  so  we 
start  by  introducing  an  expression  for  the  laminar  sublayer  thickness  -  <5  (for 

Xj 

details  see  appendix  B)  : 


1  n 
n+1  ^n+1 


-1/2 

Q 


(41) 


[appropriate  expressions  for  O  and  5  are  developed  in  appendix  A.  ] 

For  a  hydraulically  smooth  surface  G  =  12.2408  (see  the  foregoing  discussion)  and 
from  table  2,  a  =  4.63°.  Inserting  these  values  as  well  as  n  =  X  [see  Martin  (1975), 
Glass  &  Hall  (1959)]  into  equation  (39)  yields, 

xchar  =  1.25  x  10-6 


The  moan  free  path  of  the  flow  behind  the  reflected  shock  wave  at  tile  conditions 
of  the  problem  at  nand  is  10  ^  cm.  Thus  it  is  seen  that  tlie  characteristic  length 
is  of  the  order  of  the  mean  free  path. 

The  mathematical  definition  oi  a  hydraulically  smootn  surface  can  be  easily 
obtained  using  equation  (41).  At  x  =  2  cm,  which  is  half  the  length  of  our 
experimental  models: 
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As  mentioned  earlier,  a  flow  is  considered  as  hydraulically  smooth  when  the 
roughness  heights  is  smaller  than  the  laminar  sublayer  thickness.  Thus,  for  a 
hydraulically  smooth  surface 


<  6t  =  0.00517  an 
max  L 

It  is  also  of  interest  to  note  that  from  equation  (39) 
<5*  |x=2an  =  0.00086  an 


(42) 


The  obtained  values  for  xcjiar  and  emax  suggests  two  possible  criteria  for 
relating  the  roughness  characterizer  -  G  to  the  roughness  height  -  e. 

Criterion  A 


Hypothesis: 


The  above  calculated  value  of  x  ,  is  independent  of  the 

char 


roughness  height  -  e. 

Verification:  Using  the  constant  value  of  Xchar  the  values  of  G  as  a  function 
of  a  (which  depends  on  e,  table  2)  were  calculated.  Using 
this,  procedure  it  was  clear  that  for  a  given  value  of  e  they  are 
differentTBlues  of  G.  This  contradicts  our  basic  requirement 
that  C  is  independent  of  e  . 

Xchar  depends  both  on  the  incident  shock  wave  Macn  number  and 
the  roughness  height.  Therefore,  this  criterion  is  incorrect. 


Conclusion: 


Criterion  B 
Hypothesis: 

Verification: 


•me  value  off  I  2m 


height. 


6*, 


The  value  of  —  _ 

e  | x=2cm 

can  be  easily  obtained. 


6* 


x=2an 


and 


max 


is  independent  of  the  roughness 

for  a  hydraulically  smoot n  surface 
Using  the  above  calculated  values  of 
,  on.  obtains: 


Jl. 

e  x*2an  ■  0.166 
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(this  is  the  minimum  value  since  6*  was  divided  by  e  ) . 

max 

Assumption:  The  flow  over  a  rough  surface  having 

e  =  0.01  cm  belongs  to  the  frictional 

transition  flow  zone.  As  mentioned 

earlier  (equation  20) 

G  =  2.5  (InG  -  lng) 


ana  for  the  frictional  transition  flow  zone 
(equation  14)  -2.2  £  ln8  <;  0.8.  Using  an  average 
value  of  ln£5  =  0.7,  results  in: 

G  =  6.3848 

Inserting  this  value  of  G  together  with  n  =  -jl 

into  equation  (39)  results  in: 

—  _  =  0.176 

e  | x=2cxn 

Comparing  this  value  (0.176)  with  the  value 
obtained  previously  for  a  hydraulically  smooth 
wedge  (0.166)  and  recalling  that  the  latter 
value  is  the  minimum  value  for  a  hydraulically 
smooth  wedge  (these  two  values  would  be  exactly 
equal  if  e  =  0.0049  cm  which  is  smaller  than  e 


max 


was  used  for  a  hydraul ically  smooth  surface) , 


6* 

leads  to  the  conclusion  that  — 

e 

independent  of  the  rougnness  height 

Furthermore,  equation  (41)  yields: 

<5T I ~  =  0.0018  cm 

L] x=2cm 

and  since  e  =  0.01  cm: 


is 


j-  =  5.55 

°L  x=2cm 


-22- 


Conclusion: 


Summary : 


This  value  tails  well  within  the  range  1  <  £  15 

6L 

which  defines  the  frictional  transition  flow 
zone.  Thus,  the  foregoing  assumption  is 
verified . 

6* 

Let  us  now  verify  the  argument  that  —  =  const  m  the  range 

1  £  M  £  2  as  well.  For  a  hydraulically  smooth  surface  and  M 
s  s 

one  obtains: 

6*1  -  =  0.0059  cm 

x=2cm 


Dividing  this  value  by  e  =  0.0049  cm 
a  hydraulically  wedge) ,  results  in: 

=  1.204 


6*, 


e  | x=2cm 

Similarly,  for  a  rough  surface  with 


(which  is  typical  to 


e  *  0.01  an  /  one  obtains 


1.5 


ill 

e  x=2an 


1.224 


6*  i 

The  obtained  values  for  — I  „ _ 

e  I x=2cm 

is  practically  independent  of  e. 

6t | x=2cm 

JLi 


verify  our  hypothesis  that 
Furthermore, 

1.25 


—I 

e  I x=2cm 


This  value  again  verifies  our  earlier  assumption  that  the  flow 
over  a  plate  with  e  =  0.01cm  belongs  to  tne  frictional 
transition  flov.'  zone. 

6* 

The  hypothesis  that  —  is  independent  of  the  roughness 
height  -  e  is  correct. 

For  a  given  incident  shock  w dve  Kach  number 

6*i 

—  n _  =  constant 

e  | x=2cm 

In  the  framework  of  the  present  study: 


—  =  1 .224  for  1  <  M  <  2 

c  x=2cm  s 


(42) 
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6*i 

c  x=2cm 


=  0.176  for 


3.5 


£  M 


£  4.0 


(43) 


Tne  velocity  exponent  -  n 

Using  the  foregoing  criterion,  the  value  of  the  velocity  exponent  -  n  (equation 

23)  can  be  calculated  for  any  given  incident  shock  wave  Mach  number. 

Martin  (1957)  showed  that  in  the  range  1  <  Mg  <:  2  the  value  of  n  =  -^  fits 

experimental  results  best.  Glass  &  Hall  (1959)  also  arrived  at  this  conclusion. 

Based  on  these  studies  the  value  n  =  -g-  was  chosen  for  the  range  1  <  Mg  £  2. 

Schlichting  (1962)  claimed  that  the  value  of  n  decreases  when  the  flow  (shock)  Mach 

number  increases.  Trie  appropriate  value  of  n  for  any  shock  Mach  number  can  be 

calculated  using  the  fact  that  the  roughness  characterizer  is  independent  of  the  shock 

Mach  number.  For  example,  for  the  range  3.5  £  M  £  4.0  one  should  solve  tne 

s 

equation 

G(e=e. ,  M  =3.75,  n=?,  —  I  _  =  0.176)  = 

I  s  e  |  x=2an 

G(g=e.  ,  M  =1.5,  n=l  —I  -  =  1  .224). 

1  s  5'  e  |x=2cm 

Tne  solution  of  this  equation  results  in  n  =  0.1935,  in  agreement  with  the 
requirement  that  for  Mg  >  2  ,  n  <  G.2.  Tie  change  in  the  obtained  value-  of  n  is  not 
too  large.  Nevertheless,  the  present  procedure  provides  a  method  for  calculating  the 
value  of  n  for  any  incident  Mach  number.  It  is  expected  that  when  the  incident  shock 
•wave  Mach  number  increases  to  higher  values,  n  will  be  further  reduced  below  0.2. 

The  roughness  characterizer  -  C 

Using  the  values  Fr  given  in  table  2  for  any  volut-  of  roughness  height,  together 
with  the  appropriate  criterion  [equations  (42)  or  (43)]  and  equation  (39)  ,  onanleei  us 
to  find  the  values  of  the  roughness  cliaracter izer  -  C  for  a  given  roughness  height  -  e* 
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Tht  obtained  values  are  given  in  table  3.  7.s  can  be  .seen,  there  are  very  slight 

differences  between  the  obtained  value  of  G  for  a  specified  value  cf  e  .  'inis  fact 

complies  witn  our  earlier  requirement  that  the  roughness  character izer  snoufo  be 

independent  of  the  rougnness  height.  The  small  discrepancy  is  probably  due  to  tnc 

fact  thac  the  two  criteria  [equations  (42)  and  (43)]  were  established  for  Mg  =  1.5 

and  M  =3.75,  while  the  values  in  table  3  are  for  M  =  1.495  ±  0.015  and  M  =  3.74  + 
s  s  s 

0.15  • 

Tne  average  values  of  C  (given  in  table  3)  were  fitted  against  e  to  result  in 


(44) 


(  e  in  cm) .  This  curve,  as  well  as  the  (G-e)  values  given  in  table  3,  are  shown  in 
figure  15.  It  is  worthwhile  to  note  that  this  expression  complies  with  the  physical 
requirement  that  the  value  of  the  shear  stress  -t  (equation  A. 14)  increases  as  the 
roughness  height  e  increses.  Note  also  that  for  a  perfectly  smooth  surface  (e=0) , 
equation  (44)  results  in  G  ■+  °°  and  from  equation  (A. 14)  one  gets  for  this  case 


(  G  °°  )  tw  +  0  as  expected. 

A  Furtner  Cneck  of  the  Displacement  Thickness  Node! 

The  boundary  layer  displacement  thickness  model  was  based  on  experi  ••ntal  data 
obtained  for  e  =  0.01,  0.02,  0.08  and  0.2  cm.  In  order  to  consol  mate  the  findings 
additional  experiments  with  e  =  C.04  err,  were  performed. 

The  RR  ^  MR  transition  lines  obtained  by  the  present  model  for  the  ranges  1  <  Mg  £  2 
and  3.5  $  M  $  4.0  togetlier  with  the  cxjier  i.mont  ally  measured  transition  data  arc- 
plotted  on  figure  16.  The  von  Neumann  "detachment"  transition  line  is  also  snown  m 


figure  16.  It  is  clearly  seen  that  the  predict  mrs  of  tnc  present  model  are  much  letter 


than  tnat  of  tne  von  Neumann  "detachment"  criterion.  As  ixaxtai,  in  the  ram  1  f  H  £  2 

s 

excellent  agreement  is  obtained  at  the  tvnUr  of  this  range  (M  =  1.5).  Tne 

s 

predictions  of  the  present  model  for  tlie  range  3.5  i  Mg  J  4.0  are  not  as  good  as 
tho®  obtained  for  the  range  1  <  M  <  2  .  However,  tliey  arc  far  better  than  the 

O 


predictions  of  the  "detachment"  criterion. 

Conclusions 

Two  models  for  predicting  the  effect  of  the  surface  roughness  on  the  RP  t  MR 
transition  in  the  (  Mg,0w  ) -plane  were  developed.  They  are: 

1)  the  pressure  loss  model,  and 

2)  the  boundary  layer  displacement  thickness  model. 

These  two  models  suggest  alternative  ways  to  account  for  viscous  effects  when  treating 
the  RR  t  MR  transition. 

In  order  to  develop  the  boundary  layer  displacement  thickness  model,  an 
expression  for  the  boundary  layer  thickness  of  a  compressible  flow  over  rough,  moving 
flat  plates  was  developed.  Using  this  expression,  appropriate  expressions  for  both 
the  boundary  layer  displacement  thickness  and  the  laminar  sublayer  thickness  were 
derived. 

While  developing  the  boundary  layer  displacement  thickness  the  following  facts 
were  establ ished . 

1)  A  hydraulically  smooth  surface  is  a  surface  for  which  the  roughness  height 
e  <  0.00517  an 

v  <$t 

2)  The  roughness  character lzer  G  =  *  -  is  i "ceoendent  of  the  inci  dent 

v 

shock  wave  Mach  number . 

6*i 

3)  The  value  of  ~  x-2cm  does  not  depend  on  the  roughness  height. 

4)  The  characteristic  distance  -  x  ,  increases  with  decreasing 
incident  shock  wave  Pacn  number,  and  increasing  roughness  ncignt. 

5)  In  addition,  a  method  by  which  tne  velocity  exponent  -  n  can  be  determined 
for  any  incident  shock  wave  Mach  number  was  presented . 

both  models  reproduce  the  exper imental  results  guite  well. 
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Appendix  A  -  The  Boundary  Layer  Displacement  Thickness  6* 

The  origin  of  the  coordinate  system  is  moving  with  the  reflection  point  G  along 

the  wedcre  at  a  constant  ve  loci  tv  v  =  u  /cose  (u  is  the  incident  shock  wave  velocitv  in  th 

-  o  s'  w  s  - 

laboratory  frame  of  reference  and  e  is  the  wedge  angle,  see  figure  12b) .  In  this 

w 

frame  of  reference  tne  flow  in  state  (0)  is  moving  towards  the  incident  shock  wave 
(wnich  is  now  stationary)  with  a  velocity  vQ  (the  angle  of  incidence  between  this  flow 
and  tne  shock  wave  is  $  =  90  -  0^  ).  The  wall,  which  in  a  laboratory  frame  of 

reference  is  at  rest,  is  also  moving  now  with  the  velocity  =  vq  .  The  flow  in 

state  (2),  which  initially  had  an  induced  velocity  t^  (y)  ,  is  now  moving  with 

tne  velocity  V2 (y)  =  [ug  -  U2(y)]/cos0w  .  Therefore,  in  the  considered  frame  of 
reference  6*  should  be  calculated  for  a  flow  velocity  v£  (y)  over  a  moving  inclined 
flat  plate  having  a  velocity  in  the  same  direction  of  V2 (y)  * 

lor  turbulent  incompressible  boundary  layer  flows  two  different  approaches  are 
frequently  used  (both  are  semi-empirical):  they  are  the  mixing  length  of  Prandtl  and  the 
von-Karman  momentum  integrals  (see  Shames,  1982,  pp.  368).  Since  similar 
semi-empir ical  approaches  are  not  available  for  compressible,  turbulent  boundary  layer 
flows,  it  was  suggested  Dy  Eckert  (1954)  anil  Mircles  (1956)  to  adopt  the  approaches 
proposed  'hr  incompressible  flows  while  adjusting  them  to  account  for  compressibility 
effects.  In  tne  following  a  similar  approach  is  adopted. 

Prior  to  adjusting  those  approaches  to  the  compressible  case  the  differences 
between  the  incompressible  and  the  compressible  cases  should  be  understood.  While  in 
trie  former  the  velocity  and  thermal  boundary  layers  can  lx.  considered  separately,  in 
tne  hitter  tney  are  coupled.  This  coupling  arises  from  tne  fact  that  at  high  flow 
velocities  the  heat  generated  by  friction  and  the  temperature  changes  due  to 
conpre-ssibil  i ty  must  be  accounted  for.  Consequently,  in  the  compressible  case  tiic 
following  fuctors  must  be  considered; 

1)  tne  i-lacri  numlxjr , 

2)  the  Prandtl  number, 


-27- 


3)  the  viscosity  dependence  on  temperature,  and 

4)  the  heat  transfer  from  the  flow1  to  the  solid  surface. 

The  major  assumptions  upon  which  the  present  model  is  developed  are: 

1)  the  gas  behaves  as  a  perfect  gas, 

2)  the  flow  is  pseudo-steady  and  two-dimensional, 

3)  the  boundary  layer  is  turbulent  from  x  =  0.  This  assumption  is  reasonable 
because, 

a)  the  Re  number  is  very  high.  Therefore,  [Martin  (1957)]  the 
boundary  layer  is  laminar  for  a  short  distance  only,  and 

b)  Schlichting  (1962)  showed  that  the  surface  roughness 
decreases  the  value  of  the  critical  Re  number  where  transition 
takes  place, 

4)  Uie  gravitational  forces  can  be  neglected, 

3P 

5)  the  pressure  is  constant  throughout  the  entire  field,  i.e.,  —  =  0 

3P 

and  - —  =  0,  and 

3y  ' 

6)  Blasius’  semi-empirical  results  for  the  wall  shear  stress  -  r  in  inconpressible 
flews  are  applicable  when  the  average  temperature  -  Tm  is  used. 

The  v.-ill  shear  stress  -  Tw 

In  t.K-  following  an  expression  for  the  wall  shear  stress  in  a  compressible 
viscous  flow  over  a  moving  rough  flat  plate  is  developed.  The  shear  stress  for  an 
incompressible  viscous  flow  over  a  rough  stationary  flat  plate  is  (Mazor  1983 ) : 


T  2 (n-1 ) 

w  2  n+1 

—  =  u_  G 

p_ 


-  2n 


n+1 


(A.l) 


wnerc  Um  and  are  the  flow  velocity  and  density  outside  the  boundary  layer,  G  is  the 

v*6L 

roughness  cnaructerizer  G  =  — - —  (V*  -  shear  velocity,  6^  -  the  laminar  sublayer 


thickness,  v  -  kinematic  viscosity) ,  and 


^6 


u  6 
00 


V 
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Insertinq  n  =  -j  and  C  =  12.2468  (tine  value  calculated  previously  for  hydraulically 
smooth  surfaces)  in  equation  (A.1)  yields: 

—  =  0.0233  u  Rs  "1/4 

P*,  “  (A. 2) 

which  is  only  3.4%  larger  than  Blasius'  expression  for  incompressible  flows  over  a 
smooth  plate,  i.e., 

—  =  0.0225  u  Re  “1/4 
0  00 
oo 

In  order  to  apply  equation  (A. 2)  to  a  compressible  flow,  the  flow  properties  must 

be  calculated  for  the  average  temperature  inside  the  boundary  layer  -  T  . 

m 

Consequently, 


2  (n-1 ) 


—  =  u  2  G  n+1 
3m 


2n 

/%  .n+1 

u  6 


(A.  3) 


A  similar  approach  was  adopted  by  Tucker  (1951),  Eckert  (1954),  Bartz  (1955)  and 

Mireles  (1956).  By  definition:  v  =  p  /p  and  v  =  u  /p  ,  furthermore 

m  mm  °°  03  °° 

from  tlie  equation  of  state  for  a  perfect  qas:  p  =  p  /  (RT  )  , 

m  m  m 


e  =  VJWJ 


and  from  assumption  5, P  =  P 

rtn  t 


m 


Inserting  these  expressions 


into  equation  (A. 3)  results  in: 

2  (n-1) 


U!  =  u„2  G  Q 

p  00  U  6 


2n 

n+1 


(A.  4) 


wn.  v’^ 


Q  = 


2n 

,  V  n+1 

V 


(if) 

m 


1-n 

n+1 


(A.  5) 


Equation  (A. 4)  differs  as  expected  from  equation  (A.1  )  by  the  term.  Q  which 


accounts  for  compressibil ity.  In  order  to  calculate  Q,  appropriate  expressions  for 


and  T  are  required. 

m 

For  nonpolar  diatomic  gases  Mazor,  Den-Dor  &  Igra  (1985)  suggested. 


therefore 


0.64874 


T„  0.64874 
m  _  .  m. 

"  *rn  / 

O 

Inserting  this  relation  into  equation  (A. 5)  yields 

m  1-2.2975n 


Q  =  (jp-) 
m 


n+1 


Eckert  (1954)  showed  tnat  T  can  be  approximated  by; 


(A.  6) 


(A.7) 


T  =  0 .5 (T  -T  )  +  0.22 (T  -T  ) 

m  w  °°  r  00 


(A.  8) 


where  T  and  T  are  the  wall,  tiie  free  stream  and  the  recovery  temperature, 

W  IT 

respectively. 

The  recovery  temperature  -  T  (which  is  qual  to  the  wall  temperature  for  the 
case  when  there  is  no  heat  transfer)  can  be  calculated  from  (Scnliehting  ( 19G2 ) ] : 


T  u  2  u  2  Pr1/3 

-F.  =  1  +  (_w.  D  _= - 

®  00  2T 


(A.  9) 


wnere  u  and  u  are  tnc  wall  and  tree  stream  vei  ocit  ks,  reaxctively,  Ir  is  tie 

W  00 


Prandtl  numtx-r  (calculated  at  Tf)  arx:  Cp  is  tnc  sjvcifir  lent  capacity  at  constant 
pressure  (calculate  at  T  )  . 

Mirtles  (1956)  shovAxi  that  equations  (A.b)  ano  ('..9)  car.  lx  combine.'  to  result 


in : 


u  u  2 

5.56  ~  -  0.28[1  +  (— )  ] 

n  '11  '  * 


-30- 


Inscrtiny  equation  (A.  10)  into  equation  (A. 7)  results  in; 

1-2.2975n 


Q  = 


uw 

6  —  -  1 

5.56  ^  -  0.28[1  +  (^J)2]  J 


n+1 


(A. 11) 


As  mentioned,  in  the  present  case  the  frame  of  reference  is  attached  to  the 
reflection  point  G  (see  figure  12b) ,  thus:  ^  =  Vo  '  U»  =  v2  ' 


P„  =  P' 


results  in. 


,  and  v  =  v0 
oo  2 


Inserting  these  terms  into  equations  (A. 4)  and  (A. 11) 


2  (n-1 ) 


^  =  v2  G  n+1  (——) n+1 

2 


2n 


and 


v 

6  —  -  1 

n  -  J  V2 

v  -  1  - V~ - v — 2 

5.56  0.28 [1  +  (^)  ly 


(A.  12) 


1-2.2975n 

n+1 


(A. 13) 


Equation  (A. 12)  expresses  tne  wall  shear  stress  for  a  compressible  flow  over  a  rough 
flat  plate.  The  compressibility  is  accounted  for  through  C  aixl  the  roughness  is 
accounted  for  tnrougn  G.  Additional  correction  is  requried  in  order  to  use  equations 
(A. 12)  and  (A. 13)  for  the  problem  at  hand;  i.e.,  to  account  for  the  fact  that  in  the 
present  case  the  flat  plate  is  not  stationary.  This  can  easily  be  done  by  replacing 
in  equation  (A. 12),  by  the  relative  velocity  between  the  flow  and  the  surface,  V2  ~ 
vQ.  Carrying  out  this  transformation  finally  results  in: 


T  2 (n-1) 

-S-  <v2-v)2G  n+1  Q 
P2  2  .0 


2n 


/  2  \n+1 

MV2"^o|  6' 


(A. 14) 


The  boundary  layer  thickness- 6 

Inserting  the  velocity  profile  (equation  23)  and  the  wall  shear  stress  (equation  i 


.14) 
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into  von  Harman's  ntxnentun  integral  for  the  case  of  zero  pressure  gradient,  i.e., 
f~  =  0  [see  Schlichting  (1962)] 


d  r6  r  2 
Tw  »[u  - 


(v^jujdy 


results  in: 


2(n-1)  3.2975n  _  _2n_  n+1_ 

A  -  r.  3n+1  3n+1  rL-  3n+1  r(3n+1)  (2n+1),  3n+1 

x  Vry  x  1  n  J 


(A.  15) 


where  Re  = 
x 


|v2-vo|x 


Equation  (A. 15)  describes  the  boundary  layer  thickness  for  a  compressible  flow  over  a 
rough  flat  plate. 

T 

In  the  special  case  of  an  incarpressible  flow  ^  =  1,  over  a  smooth  plate 

1  A2 
(C  =  12.2468)  when  n  =  y  equation  (A. 15)  reduces  to 

,  1 

-  =  0.3816  Ite  "  5 
x  x 


Tne  expression  found  in  the  literature  for  this  case  is  ([see  Schlichting  (1962)]: 


-  =  0.376  Ve  ~ 
x  x 


2 

5 


The  difference  between  these  two  values  is  negligibly  small  (about  1.5%). 
The  boundary  layer  displacement  thickness  -{* 

The  boundary  layer  displacement  thickness  -  6*  can  bo  calculated  from 

[Schlicntmg  (1962)] 

s*  ■  s'  <’  -  r  r>  ^ 

O  oo  CO 


Inserting  tne  velocity  profile  (equation  23)  into  the  above  equation  and  carrying  out 
tne  integration  results  in: 

T 

=  6  (n+1  ■  r]  (A-16) 

m 


✓  / 


-32- 


Insorting  equation  (A. 15)  into  equation  (A. 16)  finally  yields: 


&* 

-  =  OG 
x 


2  (n-1 ) 
3n+1 


T  • 

(ft 

2 


3.2975n 

3n+1 


fcrfn) 


2n 


3n+1 


(n+1  - 


T 

-2, 
T  ’ 
m 


(A. 17) 


where 


n+1 

C  =  -L  r  <3n+1)  (2n+ll3n+l 

n+1  1  n  J  (A. 18) 

Equation  (A. 17)  describes  the  boundary  layer  displacement  thickness  for  a 
compressible  flow  over  a  rough  surface  of  a  flat  plate  which  moves  with  the  velocity v 
relative  to  the  free  stream  velocity 
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Appendix  B  -  The  Laminar  Sublayer  Thickness  -  6 

L 

The  roughness  characterizer  G  was  defined  (equation  19)  as: 


v*6l 

G  =  v 


(B.  1 ) 


where  the  shear  velocity  is: 


v*  =  (t  /p ) 


i 


w 


(B.2) 


Inserting  t  (equation  A. 14)  into  equation  (B.2)  results  in: 
w 


n-1 


v*  =  (v?-vn)Gn+1  ( 

2  0  5  v_-v, 


2  f+1  Q* 


2  ’O' 


(B.3) 


Using  this  value  (equation  B.3)  in  equation  (B. 1)  and  rearranging,  yields: 


2  1  n 

6.  -  G  n+1  (-, — - — r  n+1  <$n+1  Q_i 

L  TW 


(B.4) 


<5  and  Q  are  defined  in  appendix  A  by  equations  (A. 15)  and  (A. 13),  respectively. 
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Figure  1:  Illustration  of  four  possible  shock  wave  reflections. 

(Inter ferog rams  are  on  the  left  and  explanatory  diagrams  on  the 
right) : 

a)  regular  reflection  (RR)  ^  =  60°,  M  =  4.68. 

b)  single-Mach  reflection  (SMR)  ew  =  10°,  Kg  =  4.72. 

c)  complex-Mach  reflection  (CMR)  6W  =  20°,  Ms  =  6.90. 

d)  double-Mach  reflection  (DMR)  0W  =  40°,  Mg  =  3.76. 

The  test  gas  is  nitrogen  at  PQ  -  15  torr  and  T0  -  300K. 

( 1 , 1 1  -  incident  shock  waves,  R,R^  -  reflected  shock  wave, 

M,  M1  -  Mach  stems,  S,S^  -  slipstreams,  T,T^  -  triple  points, 

K-Kink  X  •  xi  "  triple  point  trajectory  angles, 

(0) — (5)  -  thermodynamic  states.  The  interferograms  were  taken  in 
the  UTIAS  10  x  18  cm  Hypervelocity  Shock  Tube. 

Figure  2*  Wave  configuration  of  a  regular  reflection. 

I  -  incident  shock  wave,  R  -  reflected  shock  wave, 

P  -  reflection  point,  Us  -  incident  shock  wave  velocity, 

<(>  -  angle  of  incidence,  0-  deflection  angle,  6w  -  wedge  angle. 

Figure  3:  Comparison  between  the  von  Neuman  "detachmen"  criterion  and  the 
experimental  results  on  RR  J  MR  transition  over  smooth  wedges 
in  pseudo-steady  flows.  Hie  disagreement  is  known  as  the 
von  Neuman  paradox. 

Figure  4:  Experimental  results  of  the  RR  £  MR  transition  over  rough  wedges, 
e  -  roughness  height,  and  the  "detachment"  criterion. 

Figure  5:  a)  Schematic  illustration  of  the  type  of  roughness  used  in  the 
experimental  study. 

b)  The  test  section  and  the  adjustable  test  model. 

Figure  6:  a)  SMR  over  an  e  =  0.02  cm  wedge. 

b)  RR  over  an  e  =  0.02  cm  wedge. 

c)  SMR  over  an  e  =  0.2  cm  wedge. 

d)  FP  over  an  e  =  0.2  cm  wedge. 

Figure  7:  Illustration  of  the  method  used  to  obtain  the  RR  ^  MP  transition 
wedge  angle  (8  [C  -  inverse  pressure  ratio  (Pq/P-j)  across 

the  shock  wave)! 

Figure  8:  Schematic  illustration  of  the  boundary  layer  structure  over  a 
flat  plate. 

Figure  9:  The  predictions  of  the  pressure  loss  model  against  the  experimental 
result  for  various  sizes  of  rouohness. 


Figure  10:  The  dependence  of  F  on  c. 

e 

Figure  11:  Illustration  of  the  boundary  layer  displacement  thickness  -  6*  in 
steady  flows  and  pseudo-steady  flows. 
steady  flow  :  (6*  is  positive) 

a)  velocity  profile  of  the  viscous  flow 

b)  the  equivalent  inviscid  velocity  profile. 

pseudo-steady  flow  :  (6*  is  negative) 

c)  the  unstedy  viscous  velocity  profile 

d)  the  pseudo-steady  viscous  velocity  profile 

e)  the  equivalent  inviscid  velocity  profile 

Figure  12:  Schematic  illustration  of  a  viscous  regular  reflection  at  (a) 
and  the  equivalent  inviscid  regular  reflection  at  (b) . 

Figure  13:  The  predictions  of  the  boundary  layer  displacement  thickness  model 
against  the  experimental  results  for  various  sizes  of  roughness. 

Figure  14 :  The  dependence  of  F  on  e. 

o 

Figure  15:  The  dependence  of  the  roughness  characterizer  C.  on  t  . 

Figure  16:  The  detachment  criterion  of  von  Neuman  and  the  predictions  of  the 
boundary  layer  displacement  thickness  model  for  e  =  0.04  cm  (in 
the  two  Mach  number  ranges)  and  experimental  results  of  the  RR  t  MR 
transition  for  e  =  0.04  cm. 
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PART  2 

The  reflection  of  a  planar  shock  wave 
over  a  liquid  wedge. 
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Introduction 

In  order  to  have  a  better  understanding  of  the  blast  reflection 
phenomenon  over  liquid  surfaces  such  as  seas  or  oceans  an  experimental 
study  of  the  reflection  of  a  planar  shock  wave  over  a  water  wedge  was 
carried  out. 

Experiments 

In  order  to  get  a  planar  shock  wave  reflecting  over  a  water  wedge  a 
special  shock  tube  has  been  designed  and  constructed  In  such  a  way  that 
the  whole  tube  could  be  tilted  in  a  vertical  plane.  Using  this 
technique.  It  was  possible  to  adjust  the  shock  tube  inclination  to 
obtain  any  desired  water  wedge  angle.  A  schematic  drawing  of  the  tilted 
shock  tube  with  the  water  wedge  having  an  angle  of  ew  as  well  as  a 
planar  shock  wave  reflecting  from  It  as  a  Mach  reflection  Is  shown  In 
figure  1. 

A  photograph  of  part  of  the  inclined  shock  tube  with  the  test 
section  Is  shown  in  figure  2.  A  detailed  view  of  the  test  section 
filled  In  with  water  Is  shown  In  figure  3. 

The  test  gas  throughout  the  experimental  study  was  dry  air.  Both 
the  test  gas  and  the  water  were  Initially  at  room  temperature. 

Results 

Two  different  Incident  shock  wave  Mach  numbers  were  used  with  a 
variety  of  different  wedge  angles.  The  two  Mach  numbers  were 
M1  -  1.475  ±  0.015  and  Mi  *  2.’25“±  0.05. 
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Mi  =  1.475  ±  0.015 

Typical  results  for  this  case  are  shown  In  figures  4(a)  to  4(d). 

At  ew  =  44.3°  the  Incident  shock  wave  reflects  over  the  water  wedge 
as  a  regular  reflection  [figure  4(a)].  As  the  wedge  angle  Is  decreased 
regular  reflection  becomes  Impossible  and  the  Incident  shock  wave 
reflects  over  the  water  wedge  as  a  single  Mach  reflection.  Such  a 
reflection  Is  shown  In  figure  4(b)  for  =  25.8°. 

Further  decrease  In  ew  causes  higher  triple  point  trajectory 
angles.  The  reflection  over  a  wedge  with  =  18.2°  Is  shown  In  figure 
4(c).  Figure  4(d)  Illustrates  the  case  of  a  glancing  Incidence 
<ew»o). 

It  Is  Interesting  to  note  that  In  all  the  above  photographs  a 
disturbance  Is  seen  to  be  propagating  Inside  the  water  tank.  Its 
front  propagates  along  the  air/water  interface  faster  than  the  point  of 
reflection  of  the  regular  reflection  or  the  foot  of  the  Mach  stem. 

M1  =  2.25  ±  0.05 

A  typical  regular  reflection  is  shown  in  figure  5(a)  for 
ew  =  50.8°.  At  lower  wedge  angles  0W  =  49.5°  a  double  Mach  reflection 
Is  obtained  [figure  5(b)].  At  even  lower  wedge  angles  ew  =  25.8°  the 
Incident  shock  wave  reflects  as  a  single  Mach  reflection  [figure  5(c)]. 
The  disturbance  propagating  Into  the  water  is  clearly  seen  in  this  set 
of  experiments  too. 
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The  RR  t  MR  transition  wedge  angle 

The  RR  «-  MR  transition  wedge  angles,  0 r.  ,  were  measured  for  both 

W 

tr 

M^  =  1.47  and  M^  =  2.2.  The  value  of  6W  was  obtained  by  plotting  the 

triple  point  trajectory  angle  x  as  a  function  of  the  wedge  angle  ew  for 

a  given  value  of  Mach  number.  Then  the  experimental  results  were 

tr 

extrapolated  to  get  ew  at  the  point  where  x  =  0.  The  present  results 
are  shown  In  figure  6. 

The  results  are  also  shown  In  table  1.  For  the  weaker  Incident 
shock  wave  the  actual  transition  wedge  angle  Is  about  5°  smaller  than 
that  predicted  by  the  “detachment"  criterion.  For  the  stronger  shock 
wave  the  agreement  with  the  "detachment"  transition  wedge  angle  Is 
surprisingly  good.  The  agreement  Is  better  than  that  obtained  over 
solid  wedges. 

At  the  present  a  detailed  study  of  the  reflection  over  a  water 
wedge  In  the  range  1  <  M^  <  4  Is  being  carried  out. 
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Conclusion 

The  reflection  of  a  planar  shock  wave  over  a  water  wedge  has  been 
Investigated  using  a  tilted  shock  tube  which  was  especially  designed  and 
constructed  for  this  study. 

The  four  well  known  types  of  reflection,  namely  RR,  SMR,  CMR  and 
DMR,  have  all  been  observed  over  the  water  wedges. 

The  RR  t  MR  transition  was  Investigated.  Our  early  results 
indicate  about  a  5°  disagreement  with  the  "detachment"  criterion  at 
=  1.47  and  a  surprisingly  excellent  agreement  at  M^  *  2.2.  Further 
experiments  are  being  carried  out. 


Shock  Wave 

Mach  Number 

Actual  Transition 

Wedge  Angle 

"Detachment" 

Wedge  Angle 

1.47 

43°  ±  0.5° 

48.7° 

2.2 


50.2°  ±  0.5° 


50.7° 


APPENDIX  -  LIST  OF  EXPERIMENTS  OVER  A  WATER  WEDGE 
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51.5 

2.26 

50.8 

2.30 

50.5 

2.29 

50.0 

2.31 

49.5 

2.30 

49.0 

2.29 

46.2 

2.29 

43.5 

2.21 

34.3 

2.10 

29.8 

2.22 

25.8 

2.05 
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A  schematic  drawing  of  the  titled  shock  tube  with  the  water 
wedge  having  an  angle  of  ©w. 

A  photograph  of  a  part  of  the  tilted  shock  tube  and  the  test 
section. 

A  detailed  view  of  the  test  section  with  the  water  wedge. 
Various  reflections  over  a  water  wedge  of  a  weak  shock  wave. 

a)  regular  reflection  -  M.  =  1.47  &  0  =  44.3° 

I  w 

b)  single  Mach-reflecti on  -  M.  =  1.52  &  e  =  25.8° 

i  w 

c)  single  Mach-reflecti on  -  M.  =  1.44  A  e  =  18.2° 

i  w 

d)  glancing  incidence  -  M.  =  1.47  A  e  =0 

i  w 

Various  reflections  over  a  water  wedge  for  a  moderate  shock 
wave. 

a)  regular  reflection  -  M^  =  2.30  S  8w  =  50.8° 

b)  double-Mach  reflection  -  M.  =  2.30  &  e  =  49.5° 

1  w 

c)  single-Mach  reflection  -  M.  =  2.05  A  e  =  25.8° 

i  w 

The  triple  point  trajectory  angle  -x  vs.  the  wedge  angle  -e  . 
Note  transition  from  Mach  to  regular  reflection  is  at  x  ♦  0» 
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inclined  shock  tube 


shock  wave 


